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Abstract 

The eikonal profile function J(b) obtained from the Model of the 
Stochastic Vacuum is parametrized in a form suitable for comparison 
with experiment. The amplitude and the extended profile function 
(including imaginary and real parts) are determined directly from the 
complete pp and pp elastic scattering data at high energies. Full and 
accurate representation of the data is presented, with smooth energy 
dependence of all parameters. The changes needed in the original pro- 
file function required for description of scattering beyond the forward 
direction are described. 



PACS Numbers: 12.38.Lg, 13.85.Dz, 13.85.Lg 



To appear in Physical Review D58 



1 



1. Introduction 

The application of the Model of the Stochastic Vacuum (MSV) [|l], || |3|] to 
pp and pp scattering gives a good description of the experimental quantities 
related to the amplitude in the forward direction. For large values of the 
momentum transfer \t\ the experimental amplitude has a delicate structure, 
with a characteristic dip in the differential elastic cross section that is not 
reproduced by MSV or other dynamical model without free parameters. It 
is well known that the profile function of the eikonal formalism, related to 
the amplitude through a 2-dimensional Fourier transform, has an apparently 
structureless shape, and slight changes in the profile may cause dramatic 
changes in the amplitude. The main reason is that the dominating imaginary 
part of the amplitude passes through zero for a value of \t\ which is not large, 
while the real part, being comparatively small for low \t\, fills partially the 
dip due to the zero of the imaginary part, thus producing the remarkable 
structure of the differential cross section. This is not really a shortcoming of 
MSV, since this model, in its present formulation using the first term of an 
expansion based on the eikonal formalism, has not the purpose of describing 
the differential cross section much beyond the \t\ — limit. 

Our study of the MSV profile function suggests analytical forms which 
can be used, with appropriate changes of parameters, to describe the data 
in the full range of \t\ values. In order to identify the necessary changes, 
we first describe accurately the whole data of high energy pp (at the ISR 
energies y/s= 20-60 GeV) and pp (at ^fs= 540 and 1800 GeV) scattering 
in terms of an amplitude built with a few terms controlled by a number 
of parameters. This part of our work is similar to the treatment of the 
ISR data by Carvalho and Menon [|J, but our parameters are fewer and 
have a much more regular energy dependence (which is essential to allow for 
predictions through interpolations and extrapolations). The comparison of 
the two different parametrizations exhibits some freedom of analytical forms 
allowed by the data. 

The use of an impact picture representation for the scattering amplitude 
of pp and prmp scattering was first made by Bourrely, Soffer and Wu [^]. 
With a t dependence inspired in the proton electromagnetic form factor, us- 
ing a total of 6 adjustable parameters, they are able to reproduce the general 
features of the ISR experiments. The form of the energy dependence is based 
on general properties of quantum field theory, and has allowed good predic- 
tions por the CERN measurements at 540 GeV. The difference of approach 
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with respect to our work is that we take different forms of parametriza- 
tion (inspired on the model of the stochastic vacuum) and that we do not 
parametrize the energy dependence. We show that very accurate fitting at 
every energy can be obtained with very smooth energy dependence for the 
parameters. Our description includes the data at very high \t\ (up to 15 
GeV 2 ) valid for all ISR energies and we remark that the sign of the term 
describing this tail is reponsible for differences between pp and pp data in 
this range. In the work of Bourrely, Soffer and Wu this large \t\ tail would 
be described by the Regge background term. 

We attempt to obtain from the shapes of the eikonal functions some 
understanding of the geometric characteristics of the systems in collision, 
and thus try to visualize from the data the meaning of effective radii for 
hadrons in collision |J. 



Table I - Experimental data 

Table I shows the available data on total cross sections, slope parameter 
and ratio of the forward real to imaginary parts of the amplitude in pp scat- 
tering, which come from the CERN Intersecting Storage Rings (ISR) mea- 
surements [[0]a,b,c] at energies ranging from ^fs = 23 GeV to ^fs = 63 GeV, 
and in pp scattering at the energies = 541 — 546 GeV from the CERN 
Super Proton Synchrotron (SPS) [0d] and the y/s = 1800 GeV information 
[0e] from Fermilab. In addition, there are the data on differential cross sec- 
tions in the range y/s—23 - 63 GeV [0a],||, and for pp at \fs = 546 GeV 



and y/s = 1800 GeV [10]. We also give particular attention to the behavior 
of the amplitude for large \t\ measured ]TT[ at ^J~s = 27 GeV. The data from 
Fermilab CDF [Q at 1800 GeV giving o T =80 mb seem discrepant when we 
consider the smoothness of the energy dependence of our parametrization, 
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and also is not in agreement with the predictions of Bourrely, Soffer and 
Wu ||, and thus in the present work we concentrate on the results of the 
Fermilab E-710 experiment given in Table I. 

In the MSV calculation hadronic structures enter through transverse two- 
dimensional wave-functions in the plane perpendicular to the direction of the 
colliding hadrons, for which is made the simple ansatz 



MR) = p/ir— exp (-R 2 /S 2 H ) , (1) 

OH 

where Sjj is a parameter for the hadron size. The dimensionless scattering 
amplitude Th 1 h 2 i n the eikonal formalism is given in terms of the dimension- 
less profile function Jh 1 h 2 f° r hadron-hadron scattering by 

T HlH2 =is[(g 2 FF)a*] 2 a 2 J d 2 b exp(tq-b) J HlH% {b, S x , S 2 ) , (2) 

where the impact parameter vector b and the hadron sizes Si, S2 appear 
only in units of the correlation length a, and q is the momentum transfer 
projected on the transverse plane, in units of 1/a, so that the momentum 
transfer squared is t — — \q\ 2 /a 2 . The quantity (g 2 FF) represents the gluon 
condensate. For short, from now on we write J(b) or J(b/a) to represent 
JhxHi (b, Si, S2). The normalization of T HlH2 is such that total and differential 
cross sections are given by 

a = - Im T HlH2 , — — = — — - \T HlH2 \ . (3) 
s at lbiis z 

To write convenient expressions for the observables, we define the dimen- 
sionless moments of the profile function (as before, with b in units of the 
correlation length a) 



I k = J d 2 b b k J(b) , k = 0,1,2,... (4) 
which depend only on Si/ a, S^/a, and the Fourier-Bessel transform 

I{t) = J d 2 b J (ba^f\t\) J{b) , (5) 

where jQ{ba^J\t\) is the zeroth-order Bessel function. Then 

T Hi h 2 = is[(g 2 FF)a*] 2 a 2 I{t) . (6) 
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Since J(b) is real, a T and the slope parameter B are written 



= I ll a 2 = Ka 2 . (7) 
t=o 2 J Q 

It is important to observe that these results conveniently factorize the 
dimensionless QCD strength (g 2 FF)a 4 in the expressions for the amplitude, 
and that the correlation length a appears as the natural length scale for the 
observables and for the geometric aspects of the interaction. For forward 
scattering these aspects are concentrated on the quantities I (Si/a, S 2 /a) 
and 5*1 /a, 5*2 /a) that depend on the hadronic structures and are mainly 
determined by the values of the profile functions in the range of impact 
parameter up to about 2.5 fm. The total cross section a T measures the 
strength, while the slope B has the strength cancelled out and is only related 
to the hadron geometry. The explicit formula for the slope is 

\\d 2 bb 2 Jib) 9 1 , l2 . 2 

B = =, — a 2 = - (b 2 )a 2 , 8 

2 Jd 2 b J{b) 2 w y 1 

where it is seen as related to the average value of the square of the impact 
parameter in the collision, with J(b) as weight function. We recall that here 
b is dimensionless and that (b 2 ) depends on the hadronic sizes. 

The paper is organized as follows. In Sec. 2 we present a convenient 
parametrization of the MSV profile function in terms of which the observable 
quantities are calculated. In Sec. 3 amplitudes and profile functions are 
directly obtained from experimental data and comparison is made between 
our results and those of Carvalho and Menon [|J. In Sec. 4 we present general 
comments and conclusions. 
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2. MSV profile functions 

The analytical study 0] of the MSV profile functions for large values of 
b shows that their asymptotic behavior is of the form 



J(b) = exp(-pb) 



A x A 2 

— H + ■ ■ ■ 

b b 2 



(9) 



where Ai, A 2 , ... are functions of Si /a and S 2 /a, and the quantity p = 3ti/8 
arises from the particular form of corrrelation function that has been used. 
For small and intermediate values of b the form of J(b) resembles gaussian 
shapes, which combined with the results of MSV for large b leads us to 
suggest the form 



J(b) = [J(0) - a 3 ]e- b2/ai + a 3 e- b2/a > + a 4 i 7 (6) , 



with 



Mb) 



(10) 



(11) 



The last term fixes the correct behavior of J(b) for large b given in Eq. 
@, with y/'y 2 + b' 2 written in the place of b, so that the expression can 
be extended to b=0. The combination inside the parentheses is such that 
Ay(0) = 0, and it is shown below that the use of the two terms creates a zero 
in the imaginary part of the amplitude, which is required by the data. 

The parameters a 1; a 2 , a 3 , a 4 , 7 depend on Si and 5*2. In pp and pp scat- 
tering we have Si = S 2 = S, and the behavior of the parameters as functions 
of S/a is regular and simple. Good representations are 



J(0) 



0.0024(S/a) 



4.5 + (S/a) 3 ' 

ai = 10 + 0.285 (S/af 
a 2 = 1 + 0.25 (S/a) 2 , 
0.002745(S/a) 3 



a 1 



<*3 



2187tt 2 



(S/a) 



1 + e 



( P 5/2a) 2 



(S/a) 



(12) 

(13) 
(14) 

(15) 
(16) 
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and 

7 = 0.57 {S/af . (17) 

The expression for shows the fast increase of the thickness of the asymp- 
totic tail as SJ a grows. 



FIGURE 1 

Fig. 1 - Profile function J{b) of the Model of the Stochastic Vacuum for 
S/a=2.7. The squares represent the exact values and the solid line represents 
the parametrized Jib) given by Eq. ( |l0| ) with partial contributions called Gi, G2 
and A. 

The contributions of the terms in Eq. (TTOD are illustrated in Fig. 1, 
together with the profile J(b) obtained by direct calculation [0. 

With the analytical expression (|1(]) for J(b), the integrations leading to 
the moments are easily performed, leading to 

I = 2vrQ[J(0) - a 3 K + ia 2 a 3 + , (18) 

and 

J 2 = 27rQ[J(0)-a 3 ]a? + ^a3 + ^e-^[7 + 67p]) . (19) 

As shown before |T|, 0, these expressions, which are important for the eval- 
uation of the observables, can be conveniently parametrized as functions of 
S/ a, leading to the forms 

J = 0.006532 (S/a) 2 - 791 (20) 
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and 



1 h 

K = ~ = 2.030 + 0.3293 (S/a 

2 Jo 



,2.126 



(21) 



The evaluation of the amplitude through Eqs. (0) and (|6|) requires the inte- 
gration formula 



W U )—f 2 , | 



- 7v /p 2 +/3 2 



Vp 2 + P 2 



and then we obtain 

I(t) = 2vr 
where 

A,(t) 



l[J(Q) - a 3 } ai e~ a2 ^ + ^a 3 a 2 e 



-a 2 \t\a 2 /i 



(22) 



(23) 



Ay(b)J (baJ\t\)bdb 



P 2 +a 2 \t\ 

p 2 + a 2 \t\ 



— e 



-l^J^P 2 +a 2 \t\ 

4p 2 + a 2 \t\ 



(24) 



The difference of two terms in this expression accounts for a zero in the 
imaginary amplitude. With the physical value a = 0.32 fm, the dimensionless 
combination a 2 |t| takes the reference value a 2 \t\ = 1 for \t\ = 0.38 GeV 2 . For 
large \t\ ( for a 2 \t\ ^> 1) the largest contributions to the amplitude I(t) come 
from the last term, with 



I(t) -> 2tt a 4 



--Y\/a 2 \t\ 



:i-e 



IP 



2, a 2 \t\ 



1 - 4e 7p ) + . . . 



(25) 



and we observe that for large \t\ the amplitude is negative (7/9 is larger than 
1). For small \t\ the behavior is 



7(t) «27r|i[J(0)-a 3 ]aie- 



« 2 l*l«i/4 + £ -a*|i|a 2 /4 + a i 

2 3 2p 



1 



37 7 
4> + 8p^ 



cr t + ... 



(26) 



where the last two terms are small compared to the first, so that the am- 
plitude has the typical exponential behavior. However, we remark that the 
amplitude is not of the simple form T(s, \t\) = T(s, 0)e _B '''/ 2 , and the slope 
B must be evaluated through Eqs. (|7|), (||), ([19]). For t=0 we have I = 1(0) 
consistently with Eq. (pT8|) . 

The results obtained with MSV are correctly related to experimental data 
at t=0 (total cross section and slope), leading to the determination of the 
QCD parameters |l|, [|. However, the predictions for the differential cross 
sections beyond the small values of \t\ are not correct, mainly because the 
MSV profiles put the zeros of the imaginary amplitudes at too large values of 
\t\. In the next section we introduce empirical modifications in the parameters 
in order to reproduce the data in the full \t\ range. 

3. Amplitudes and profile functions from the experimental data 

Using the same form of parametrization suggested by MSV, we construct 
the dimensionless amplitude 

T cxp (s, t) = A^s[iX{t) + K{t)) (27) 

directly from the experimental data (hereafter, expressions such as experi- 
mental amplitude, experimental profile and so on will be used for quantities 
obtained directly from the data). The imaginary and real parts of the am- 
plitude are parametrized using the forms 

X{t) = ai e~ m + a 2 e~ m + A 2pe P7 A y (t) (28) 

and 

K{t) = X' 2pe^'A Y (t) + a[e-^ , (29) 

where p = 3n/8 and where we have grouped factors 2pe P7 A 7 (t) in order to 
have 2pe P7 y4 7 (0) = 1. We use 8 parameters for each energy (here we do not 
count a[ and (3[ which are the same for all energies, as explained below). 

The use of Ay(t) in the real part, which is inspired in the MSV form for 
the imaginary part, has a more convenient structure, compared to simple 
exponentials, to fill the dip left by the zero of the imaginary part. On the 
other hand, the simple exponential term a[ exp(— f3[\t\) was included in the 
real part specifically to describe the large \t\ (5 < \t\ < 15 GeV 2 ) data ]TTJ| at 



27 GeV, which is assumed to be universal for all ISR energies. This term is 
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not present in our description of the data at 546 and 1800 GeV, where large 
\t\ values have not been measured. 

We use as experimental inputs for pp scattering at ISR energies and for 
pp at 546 and 1800 GeV the ratio of the forward real to imaginary parts of 
the amplitude, 1Z(0)/X(0) = p, the slope at t = 0, 



dt\ dt] 



K(0)K'(0) +T(0)T'(0) 
t =0 ft 2 (0)+Z 2 (0) ' 1 ] 



and the total cross section determined through the optical theorem, 



a 



4^J(0) . (31) 



After this information is introduced, the remaining parameters are deter- 
mined by fitting \TZ(t) + iZ(t)| 2 to the experimental data for da/dt. The 
results of fittings at the ISR energies and at 546 and 1800 GeV are shown in 
Figs. 2 and 3. 



Fig. 2 - Results of fittings of elastic differential cross sections through Eqs. 
(f2~§|) and (^9|) at ISR energies. Curves and data are conveniently separated in the 
figure through multiplication by the factors that are shown in the upper corner. 
The large \t\ data measured at 27 GeV that are drawn together with all angular 
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distributions are described by the last term of the real amplitude in Eq.(2£). 



FIGURE 3 

Fig. 3 - Same as in Fig. 2 , for 546 GeV and 1800 GeV. 

The parameters present a very regular energy dependence, as shown in 
Table II and Fig. 4. However, it is dificult to write parametrized analytical 
dependences for the parameters as functions of \/s because the expressions 
given by Eqs. (28) and (|29|) are very sensitive to the values of the parameters, 



requiring accuracy of at least three decimals to reproduce correctly the data 
if we require smoothness in the variation with the energy. If each energy is 
considered separately, and allowing for the experimental errors in the data, 
the parameters are fixed typically within 1% at the energies for which the 
data are more accurate (e.g. 52.8 GeV). 



Table II 

Table II - Values of the parameters for Eqs. (|28| ) and fl29|). 
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FIGURE 4 

Fig. 4 - Parameters of the fittings through Eqs. ( f28|) and ( |29|) at ISR energies, 
and at 546 and 1800 GeV. The energy dependence of the total cross section a T 
and the slope B obtained from the parameters are also shown. 

The imaginary and real parts of the experimental scattering amplitude 
T e xp{s,t) are shown in Figs. 5 and 6 respectively. The figures show remark- 
able regularities in the energy dependence, with an interesting similarity in 
the behavior of all curves as functions of \t\. The real parts grow with the 
energy faster than the imaginary parts, so that even after multiplication by 
1/s the real parts still have magnitudes which increase with s (note the in- 
verted ordering of the curves in the two figures). At all investigated energies 
the imaginary amplitudes present only one zero as functions of \t\, with a 
regular displacement of the zeros towards the origin as the energy increases 
(shrinking of the diffraction peak). The prediction of a second zero in the 
imaginary part by Borrely, Soffer and Wu || for very high energies ( above 
10 TeV) is very far for the range of our analysis. Depending on the sign 
given to the last term in Eq. (|29"D, the real amplitude may present one (for 
negative sign) or two zeros (case of positive sign). This last term practically 
does not interfere with the first part, as its contribution is only important in 
the tail of the distribution, and thus the fitting does not determine its sign. 
The first zeros of the real parts are distributed in a rather small range of 
\t\ values, closer to the origin compared to the zeros of the imaginary parts. 
The displacement of the zeros with the energy is shown in Fig. 7 (only the 
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first zeros of the real parts are shown). 



FIGURE 5 

Fig. 5 - Imaginary part of the experimental scattering amplitude at the ISR 
energies and at 546 and 1800 GeV, divided by s. 



FIGURE 6 

Fig. 6 - Same as in Fig. 5, but for the real part of the amplitude. 



FIGURE 7 
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Fig. 7 - Locations of the zeros of the imaginary and real parts of the amplitude 
as functions of the energy. 



Unfortunately the data for da/dt at 1800 GeV do not go beyond \t\ = 
0.6 GeV 2 , and the positions of the zeros and the large \t\ behavior cannot 
be extracted from the data. Actually, Figs. 5, 6, 7 only give (extrapolated) 
guesses for the quantities related to this energy. 

Data at large \t\ (up to 15 GeV 2 ) obtained [T| only at 27 GeV seem to 
offer a natural extrapolation for the angular distribution at all ISR energies, 
as shown in Fig. 2. This tail is described by the term c^e - ^ 1 '*' in the real 
part in Eq. (|2D|), with small values for a[ and j3[, which hardly influences 
the determination of the other parameters. Thus we have a term specific for 
the energy independent tail, which has been included in the curves of Fig. 
2. The large \t\ tail in da/dt has been parametrized [D| in the form 0.09/t 8 , 
and to have an expression with this behavior and that at the same time can 
be used for all \t\, we may modify Eq. fl29|) to have the form 

(32) 



K{t) = X' 2pe^A yl (t) + a'[[l - e -/WI*l) 4 ] 



We discuss this representation in the next section. 

To obtain the profile functions from the experimental data, we invert the 
Fourier-Bessel transform 

(33) 



n{t) + il(t) = 2tt / [£(&) + iJ(b))J (bJ \t\)bdb 



where b = ab (here b is in fermi or GeV 1 units). We then obtain the desired 
(dimensionless) imaginary and real parts of the profile function 



J(b) 



and 



1 

£(&) = 



1 

47T 



-b 2 /A(3 



1 + 



^2 -6 2 /4ft 



+ 4^e^i 7 (6/a) 



A Y (b/a) + ^e-' b2 ^ 



(34) 



(35) 



To follow the same notation used for the profile function of MSV (see 
Eq. (|10D and recall that in MSV we have only the imaginary part of the 
amplitude) we must introduce a numerical factor 



Jexp(b) 



4^ 



[{g 2 FF)a A \ 2 



Jib) 



(36) 
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with A^/[(g 2 FF)a 4 } 2 = 2.02 x 10~ 2 . 

Given the experimental profiles JC(b) and J~{b), we can now write con- 
venient expressions for comparison between the experimental quantities and 
those of MSV. We define the moments of the experimental profile functions 
corresponding to the imaginary part of the scattering amplitude, 



X fe = 2tt / b k J(b)bdb , (37) 



and to the real part 



TZ k = 2vr / b k IC(b)bdb , (38) 



Jo 



which depend only on y/s. The moments can be evaluated analytically, and 
are given by expressions similar to Eqs. ( |T8"D and (fOJ). 

Including the real part, the expression for the slope parameter at t — 
in terms of the moments TZk and is 

TZ TZ 2 +1 1 2 pTZ 2 +l 2 
^ cxp 2(72g+2g) 2(l+p2)J ' ^ 

where TZq/Xq = p . The maximum error in -B eX p caused by the suppression 
of the real part contribution is about 2% at 1800 GeV and less than 1% at 
the ISR energies. 

In order to compare experimental quantities with those of MSV, we must 
have a correspondence between the experimental values of y/s and the MSV 
values of S/ a for the proton radius |], . We then require that 

W'» =3i( ' ) [tf^v (40) 

which means that the evaluated total cross sections must be the same. The 
solutions of this equation can be written in a simple parametrized form 

S/a= 2.1848 + 0.1623 In Vs. (41) 

According to the phenomenology of the application of MSV to pp and pp 
scattering [|I|, [J, the experimental and MSV values of o T and B are fixed to 
be the same for each one of the two input energies (541, 1800 GeV). This 
fixes the two fundamental QCD quantities (gluon condensate and correlation 
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length) and the values of the effective proton radius at each of these two 
energies. Then the MSV and experimental calculations of I2 (and obviously 
the ratio I2/I0 that gives the slope B) agree within 1 to 2% accuracy at all 
energies. This is the basis of the success of the MSV calculation of total cross 
sections and slopes. 

It is important to compare our results with those obtained by P. A. Car- 
valho and M. J. Menon (CM) who parametrized the scattering amplitude 
(here rewritten according to our notation) using a sum of exponentials in the 
form 

2 n 

F(s, t)= - n(s) £ a 3 e-^ + i ]T aye"^* 1 , (42) 
j=i 3=1 

where 

with coefficients obtained by fitting da/dt = 7r\F(s,t)\ 2 to the pp experi- 
mental data at the ISR energies (n= 4-6, depending on the energy). The 
corresponding dimensionless profile functions are 

KcM(6) = -^Efe-^ (44) 
A j=i Pj 

for the real part, and 

JcM(b) = W^- b2 ^ (45) 
z 3=1 Pj 

for the imaginary part. Comparison with our formulas is given by the corre- 
spondence of values 

/Com® = 2v^/C(6) , J cu {b) = 2^J{b) . (46) 

Our imaginary amplitudes Im[T exp (s, £)]/4a/7ts = Z(i) and the val- 
ues of y/7rlm[F(s, t)] obtained by Carvalho and Menon are very similar for 
all ISR energies. The Fourier tranforms J{b) and J7cm(&), however, coin- 
cide only between b=0 and b=2. On the other hand the real amplitudes 
Re[T exp (s, t)]/4y/7rs = 7Z(t) and y / 7rRe[F(s, t)] differ substantially for all val- 
ues of \t\, except at \t\ = where they coincide (within the experimental 
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error), because of the input experimental value for p. The comparison of the 
real parts obtained in the two works is illustrated in Fig. 8 for the amplitudes 
and in Fig. 9 for the profile functions. 



FIGURE 8 

Fig. 8 - Comparison of the real parts of the experimental amplitudes at 23.5 
GeV obtained in our work (solid curve) and by Carvalho and Menon (CM) H 
(dashed curve). 



FIGURE 9 

Fig. 9 - Same as Fig. 8, but for the real parts of profile functions, with the 



correspondence between profiles given in Eq. (46). In our representation (solid 
line) the value of the profile function at b = is fully given by the second term of 
the amplitude, that only refers to the tail. 

In spite of these differences in the real parts, the two parametrizations are 
equivalent with respect to the description of the data at each separate energy. 
However, our results are distinguished by the very regular energy dependence 
of the parameters. This regularity, which is important for interpolations and 
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extrapolations, is true up to the highest energies 540 and 1800 GeV and also 
includes the angular distribution at large \t\ measured at 27 GeV. 

4. Comments and Conclusions 

The application of MSV to high-energy scattering leads to a good descrip- 
tion of the observables determined by the forward amplitude (cr T and B), but 
the predictions for the differential cross sections for large \t\ are not correct, 
mainly because the MSV profiles put the zero of the imaginary amplitude at 
values of \t\ that are systematically higher than those required by the data. 
The model gives a dynamical, QCD based, framework to evaluate the forward 
amplitude in the eikonal formalism, and in its present formulation does not 
intend to describe elastic scattering at all \t\, that would require the inclusion 
of a real part in the amplitude and changes in the imaginary part. In the 
present work, starting from the analytical forms of profile function and am- 
plitude suggested by the model, we investigate the modifications needed to 
fit all elastic scattering data of the pp and pp systems at high energies. Val- 
ues are obtained for the parameters and we believe that the results are made 
rather unique through the effort to obtain a smooth energy dependence. The 
chosen form respects the properties of the profile function and is such that 
the 2-dimensional Fourier transform can be written analytically. The exis- 
tence of a zero in the imaginary amplitude is a common feature of both the 
MSV calculation and the experimental data, and the measured differential 
cross section requires that this zero be filled up with the contribution from 
the real part. We thus must add a real part to the MSV amplitude. 

In our parametrization the 27 GeV tail (5£ \t\£ 15 GeV 2 ) is specifically 
described by a term in the real part. The large \t\ amplitude is probably dom- 
inated by perturbative contributions, such as triple gluon exchange, which 



is real and can account for the observed 0.09/t behavior [fL3|, although it is 
surprising that there is no energy dependence. Fig. 10 shows the compari- 
son between the exponential form used in Eq. (p9|) and the power behavior 



produced by Eq. (R2T). In this last case the values of the parameters are 
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a'l = 0.001 GeV" 2 and f3'{ = 23650. Both forms describe well the data. 



FIGURE 10 

Fig. 10 - Differential cross section for large \t\ values. Comparison between 
the parametrization given by Eqs. (28) and (|29|), and the l/|i| 8 dependence [13 
produced by Eq. (|3 



In Table II we have chosen a positive sign for the tail term added to 
the real part, although it could also be negative, the different choices only 
requiring a slight adjustment of parameters. A positive sign in this term 
causes a second zero in the real amplitude occurring at \t\ « 2 GeV 2 , but 
there appears no dip because the imaginary part is large in this region. The 
first zero of the real amplitude occurs at very small \t\ for either sign of the 
tail term. The roles of the real and imaginary parts of the amplitude, and 
separately of the real term for the tail, considering its two possible signs, are 
shown in Fig. 11. 



FIGURE 11 

Fig. 11 - Separate roles of |Im T| 2 and |Re T| 2 at 23.5 GeV. The dip in the 
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data at about 1.5 GeV 2 is due to a zero in Im T. The real part is positive at t = 0, 
becomes negative at small \t\ and fills partially the dip in da/dt. The dotted line 
represents the squared tail term alone. If this term is positive, it causes a second 
zero when added to Re T, as shown by the (+) line. In the (-) line the tail is 
negative and no second zero exists. 

Fig. 12 shows the profile functions obtained from MSV and from the 
fitting to the data. A general fact is that the experimental profile functions 
are lower for small b and higher for intermediate b values, compared to the 
MSV shapes. Due to the product bJ(b) that appears in the 2-dimensional 
Fourier transform, the region about b ~ has little physical significance. It is 
then more instructive to examine the forms of the product bJ(b), as shown in 
Fig. 13. We remark that in this figure corresponding MSV and experimental 
profiles have the same area (the moments Iq fix the total cross section). In a 
plot of b 3 J(b) the corresponding profiles also have the same areas, since the 
MSV moments I2 yield the correct slopes. 



FIGURE 12 

Fig. 12 - Comparison of values of J(b) from experimental data and from MSV, 
for corresponding values of ^/s and S/a according to the requirement established 
byEq. ©. 
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FIGURE 13 

Fig. 13 - Same as in Fig. 12, but for b,J(b). 



The present development of the MSV formalism applied to hadron-hadron 
scattering [Q is restricted to the lowest order nonvanishing contribution, 
which is quadratic in the gluonic correlator, and has been shown to be suf- 
ficient for the evaluation of total cross sections and slope parameters. The 
resulting amplitude is purely imaginary, and the ratio of the real to the 
imaginary parts of the elastic scatering amplitude can only be described if 
we go one further order in the contributions of the corrrelator. On the other 
hand, we have verified that changes in the form of the correlation functions 
D(z 2 /a 2 ) and Di(z 2 /a 2 ), that enters the definition of the correlator, and rea- 
sonable changes of the hadron wave-function, have little effect on the results. 

In addition to these considerations, which concern mainly the nonpertur- 
bative nature of the model, we must also take perturbative contributions in 
account. To the MSV calculation, which represents single pomeron exchange, 
we may add terms with mixed pomeron and multiple gluon exchanges, which 
modify the imaginary and real parts, and may add three gluon exchange 
(purely real) for the tail JTH . 



In Fig. 14 we plot the difference at 23.5 GeV between the experimental 
and MSV imaginary parts of the amplitude T(s, t) 

- [Im T exp (s,t) - Im Tmsv(s,£)] , 

which illustrates the contribution to be added to the MSV result for the imag- 
inary part, which must be mostly negative, peaked at about \t\ = 0.5 GeV 2 . 
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In the same figure we plot the real part that is extracted from the data. We 
observe that, except for a small range of very small \t\ values, the corrections 
to the imaginary part must be much larger than the values of the real am- 
plitude. The situation is similar for the other energies. These remarks may 
guide future theoretical developments. 



FIGURE 14 

Fig. 14 - Real part of the amplitude and difference between the imaginary 
parts of the amplitude obtained from the experiments and from the MSV calcula- 
tion, all at 23.5 GeV. 

Fig. 15 shows the real parts of the profile function for all energies of our 
investigation. We observe a very regular energy dependence. The quantity 
Ay'ib/a) in Eq. ( |35|) is zero for b=0, so that the value at the origin is due 
only to the tail term, being the same for all ISR energies, and not included 
at 540 and 1800 GeV. 



FIGURE 15 

Fig. 15 - Real part 1C(b) of the experimental profile functions for all ISR 
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energies, 546 GeV and 1800 GeV. 



Now we comment on the relation between observables in high-energy scat- 
tering and the hadronic sizes. The geometrically extended character of the 
non-perturbative QCD interaction determines the phenomenological proper- 
ties of the observables, which are fixed by the sizes and global structures of 
the colliding systems. These features have led to models of geometric nature 
for high-energy scattering |J, that give a natural account of the relations 
between total cross sections and slopes of the diffractive peaks of different 
hadronic systems. A geometrical relation of this sort is seen if we write Eq. 
(P[) in the form 

B = 1 -i^BB= 1 -m. (47) 

2 fdb[bJ(b)] 2 W K J 

and interpret (b 2 ) as the average value of the square of the impact parameter 
taking bJ(b) as weight function. There is an interesting relation between this 
average value (b 2 ) and the square of the value b$ of the impact parameter 
at the peak of the function b 3 J(b) that appears in the upper integrand of 
Eq. (f4"Tp. We have verified that the ratio [&£P] 2 / (b ) of these two quantities 
extracted from the data is a constant in a wide energy range, as shown in 
Fig. 16. In the case of a pure gaussian profile (i.e. an amplitude with a 
single exponential in t) this ratio is 1.5 . This property help us to define an 
effective hadronic radius, related to b$, in terms of a measurable quantity. 
It is remarkable, and may be of practical importance for the determination 
of slopes, that the position b$ of the peak is not dramatically sensitive to 
the details of the amplitude close to t — 0. 



FIGURE 16 
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Fig. 16 - Values of the experimental ratio [bmlx] 2 /(b 2 ) (black squares) and its 
mean value (line). 

A convenient parametrization of the experimental values of is 

fegk = 1.0119 + 0.0433 hiv^ , (48) 

and with the correspondence between S/a and y/s given by Eq. (fl~ID , we 
may also write 

ggk = 0.4290 + 0.2668 (S/a) . (49) 

The quantity fe^L (written above in fermis, while yfs is in GeV) can 
be extracted from the amplitude and gives a measurement of the effective 
hadronic size. Eq. ( fi9|) connects the effective hadronic collision radius with 
the radius of the transverse wave-function. As the energy varies from 23 to 
1800 GeV, 6gk varies from 1.15 to 1.34 fm. We obtain for yfs=U TeV the 
predictions = 1.43 fm, S/a = 3.73, and B exp = 18.76 GeV" 2 . 
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Table 1: Experimental Data 







a 1 


B 


Ref. 


P 


Ref. 




(GeV) 


(mb) 


(GeV 2 ) 


i 




01 




23.5 


39.65 ±0.22 


11.80 ±0.30 


(a) 


0.020 ± 0.050 


(b,c) 




30.7 


40.11 ±0.17 


12.20 ±0.30 


(a) 


0.042 ±0.011 


(6,c) 


pp 


44.7 


41.79 ±0.16 


12.80 ±0.20 


(6) 


0.062 ±0.011 


(b,c) 




52.8 


42.38 ±0.15 


12.87 ±0.14 


(a) 


0.078 ±0.010 


(b,c) 




62.5 


43.55 ±0.31 


13.02 ±0.27 


(a) 


0.095 ±0.011 


(b,c) 


pp 


541 


62.20 ± 1.50 


15.52 ±0.07 


(d) 


0.135 ±0.015 


(d) 




1800 


72.20 ± 2.70 


16.72 ±0.44 


(e) 


0.140 ±0.069 


(e) 



Table 2: Values of parameters for Eqs. ( |28| ) and ( p9[ ) 



>/5 






a 2 




A 


A' 


7 


y 


x 2 


23.5 


1.7298 


1.4390 


3.1091 


2.2949 


9.1850 


0.2774 


3.8620 


6.00 


1.16 


30.7 


1.8224 


1.4502 


3.1649 


2.3299 


9.5467 


0.6073 


3.9000 


6.600 


1.21 


44.7 


1.6699 


1.5025 


3.0000 


2.1086 


10.3650 


0.9354 


4.0311 


6.550 


4.28 


52.8 


1.8500 


1.5287 


2.9600 


2.1753 


10.4630 


1.1913 


4.0540 


6.685 


3.22 


62.5 


1.9272 


1.5529 


2.8081 


2.1337 


10.8201 


1.4747 


4.0924 


7.300 


1.75 


546 


2.6174 


2.1539 


3.9061 


2.1813 


15.7463 


3.0270 


4.8190 


7.970 


1.23 


1800 


3.1036 


2.4526 


4.4246 


2.4253 


18.3315 


3.6204 


5.3450 


8.600 


2.00 


a[ = 


«i, Pi, «2, P2, A and A' are in GeV 2 . 7, 7' are dimensionless. 

= 0.0031 GeV" 2 and p[ = 0.41 GeV" 2 are the same for all ISR energies. 
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FIGURE CAPTIONS 



Fig. 1 - Profile function J(b) of the Model of the Stochastic Vacuum for 
S/a=2.7. The squares represent the exact values [|J, and the solid line represents 
the parametrized J(b) given by Eq. ( |l0| ) with partial contributions called G±, G2 
and A. 

Fig. 2 - Results of fittings of elastic differential cross sections through Eqs. 
(p£|) and (^9|) at ISR energies. Curves and data are conveniently separated in the 
figure through multiplication by the factors that are shown in the upper corner. 
The large \t\ data measured at 27 GeV that are drawn together with all angular 



distributions are described by the last term of the real amplitude in Eq.(2£). 
Fig. 3 - Same as in Fig. 2 , for 546 GeV and 1800 GeV. 



Fig. 4 - Parameters of the fittings through Eqs. ( pq ) and (|29|) at ISR energies, 
and at 546 and 1800 GeV. The energy dependence of the total cross section a T 
and the slope B obtained from the parameters are also shown. 

Fig. 5 - Imaginary part of the experimental scattering amplitude at the ISR 
energies and at 546 and 1800 GeV, divided by s. 

Fig. 6 - Same as in Fig. 5, but for the real part of the amplitude. 

Fig. 7 - Locations of the zeros of the imaginary and real parts of the amplitude 
as functions of the energy. 

Fig. 8 - Comparison of the real parts of the experimental amplitudes at 23.5 
GeV obtained in our work (solid curve) and by Carvalho and Menon (CM) Q 
(dashed curve). 

Fig. 9 - Same as Fig. 8, but for the real parts of profile functions, with the 



correspondence between profiles given in Eq. (46). In our representation (solid 
line) the value of the profile function at b = is fully given by the second term of 
the amplitude, that only refers to the tail. 

Fig. 10 - Differential cross section for large \t\ values. Comparison between 
the parametrization given by Eqs. (28) and (|29|), and the l/|i| 8 dependence [13] 



produced by Eq. (p 
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Fig. 11 - Separate roles of |Im T| 2 and |Re T\ 2 at 23.5 GeV. The dip in the 
data at about 1.5 GeV 2 is due to a zero in Im T. The real part is positive at t = 0, 
becomes negative at small \t\ and fills partially the dip in da/dt. The dotted line 
represents the squared tail term alone. If this term is positive, it causes a second 
zero when added to Re T, as shown by the (+) line. In the (-) line the tail is 
negative and no second zero exists. 

Fig. 12 - Comparison of values of J(b) from experimental data and from MSV, 
for corresponding values of ^fs and S/a according to the requirement established 



by Eq. (g0|) 



Fig. 13 - Same as in Fig. 12, but for bJ(b). 

Fig. 14 - Real part of the amplitude and difference between the imaginary 
parts of the amplitude obtained from the experiments and from the MSV calcula- 
tion, all at 23.5 GeV. 

Fig. 15 - Real part fC{b) of the experimental profile functions for all ISR 
energies, 546 GeV and 1800 GeV. 

Fig. 16 - Values of the experimental ratio [&max] 2 /(fr 2 ) (black squares) and its 
mean value (line). 
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